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A superconductor with interacting paramagnetic impurities is considered. The impurities are 
coupled via the Ruderman-Kittel-Kasuya-Yoshida interaction. At a temperature T g , the system of 
magnetic impurities forms a spin-glass state. We study the effect of the spin-spin interactions on the 
superconducting transition point at T < T g . We show that superconducting properties depend on 
the state of the spin system via spin-spin autocorrelation functions. With the help of the Keldysh 
technique, a general nonequilibrium Gor'kov equation is derived. Possible ageing effects in the 
superconducting transition point are discussed. The equilibrium superconducting transition point is 
found explicitly and shown to be shifted towards higher temperatures and impurity concentrations 
compared to the classical Abrikosov-Gor'kov's curve. The corresponding shift of the superconducting 
quantum critical point is quite small (about 10%). A method of calculating spin-spin correlation 
function is suggested. The method combines the ideas of random mean-field method and virial 
expansion. We calculate analytically the first virial term for the spin-spin correlator for the quantum 
Heisenberg spin glass with the RKKY interactions in the quasiequilibrium regime. 
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I. INTRODUCTION 

The theory of superconducting alloys with paramag- 
netic impurities was developed long ago by Abrikosov 
and Gor'kov .El They have shown that the superconduct- 
ing transition temperature was suppressed by magnetic 
impurities. At some critical impurity concentration, the 
transition temperature was suppressed down to,-zero, 
which gives an example of a quantum critical point .□ The 
critical concentration can be determined from the condi- 
tion: t s T C {) = 27/71", where T c q is the superconducting 
transition temperature in a sample without impurities 
and t s is the spin-flip scattering time, which is inverse 
proportional to the concentration of magnetic impurities 
n s . Let us note that the effect of magnetic impurities on 
superconductivity is, in manv_.aspects, similar to the one 
of an external magnetic fieldH 

In the Abrikosov-Gor'kov's model, the magnetic impu- 
rities did not interact. Indeed, in any real system, one 
can not avoid having interactions between the impurities. 
Friedel oscillations in the electron density give rise to the 
similar oscillations in the spin-spin coupling which is well- 
known as Ruderman-Kittel-Kasuya-Yoshida (RKKY) in- 
teraction. The interaction may change the nature of the 
transition both qualitatively and quantitatively. 

Interactions between impurities can significantly 
change the physical picture only if the typical energy of 
this interaction is of the order of temperature or higher. 
At a high temperature, one can take into account only 
those impurity clusters in which the typical distance be- 
tween spins is small enough. The probability of having 
such a cluster containing three or more spins is small. 
Thus, virial expansion (or cluster expansion) coincides 
with the high-temperature expansion in the problem un- 



der discussions The corresponding expansion parameter 
is Jqti s /T, where Jo is the amplitude of the RKKY in- 
teraction. 

Let us note that constant Jo and spin-flip scattering 
time t s both follow from the same exchange Hamiltonian. 
Parametrically, they are of the same order. However, due 
to the numerical smallness of the two-particle statistical 
weight in the three-dimensional case, the amplitude of 
the RKKY interaction is numerically much smaller than 
the spin-flip scattering rate: JoT s n s — 1/4it 2 S (S + 1), 
where S is the spin of a magnetic impurity. 

The high-temperature corrections to the Abrikosov- 
Gor'kov's result were derived in the paper of Larkin, Mel- 
nikov and KhmelnitskiiH It was shown that the supercon- 
ducting transition temperature was higher in the pres- 
ence of the spin-spin interactions compared to the non- 
interacting case. However, the corresponding change was 
noticeable only at extremely low temperatures due to the 
small numerical factor mentioned above. Thus, in a very 
wide range of temperatures, the Abrikosov-Gor'kov's the- 
ory was quantitatively correct. Let us mention that there 
have been a number of experiments in which deviations 
from the Abrikosov-Gor'kov's picture were observed at 
low temperatures (see e.g. Ref. ||). 

Since the RKKY interaction is random in sign, it intro- 
duces a frustration into the spin system. It may result 
in a spin glass transition at a low enough temperature 
Tg ~ Jo^s- This makes the physical picture much more 
puzzled compared to the high-temperature case. A wide- 
spread distribution of energy barriers exists in the glassy 
phase. The typical time of classical and quantum tunnel- 
ing is comparable with the observation time in real exper- 
iments. Thus, the state of the system depends not only 
on the temperature and external magnetic field, but also 
on the history. Moreover, all physical quantities slowly 
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depend on time. OS 

In the present paper, we show that superconducting 
properties, transition temperature T c in particular, de- 
pend on the properties of the spin system. Thus, super- 
conducting parameters are expected to depend on the 
history and real time if the spin system is in the ageing 
regime. 

Although the dynamical properties of the spin glasses 
have been the subject of extensive studies in the recent 
years, the theoretical understanding of the effect is far 
from being impressive yet. However, it is clear that, in 
general, the dynamics of a glassy system consists of two 
parts: one is a fast quasiequilibrium dynamics and the 
other is a slow dynamics or ageing. One can expect that 
in the system under consideration different superconduct- 
ing parameters should acquire similar behavior. 

Studying the out-of-equilibrium dynamics in a quan- 
tum system requires using a nonequilibrium formal- 
ism such -as quantum transport equation or Keldysh 
technique B The Keldysh technique in superconductors 
was developed bj£-.Larkin and OvchinnikovLj and by 
Feigel'man et al!El. In the theory of quantum spin 
glasses, Keldysh formalism was recently used by Kennett 
et al£3 In the theory of classical spin glasses, Keldysh, 
technique is replaced by so-called Doi-Peliti techniques ,t2l 
in which dynamics is generated via the Langevin noise in- 
troduced in the stochastic equation of motion. This kind 
of technique was used by Ioffe et, al. who considered dy- 
namics of a classical spin glass. L3 

Our paper is structured as follows: 
In Section II, which is necessarily quite technical, we 
derive general equations on the superconducting Green 
functions taking into account a possible non-equilibrium 
dynamics in the spin system. The corresponding cal- 
culations are done with the help of the Keldysh tech- 
nique. We reconsider the Abrikosov-Gor'kov's theory, 
taking into account inelastic exchange electron scatter- 
ing on magnetic impurities. Technically speaking, impu- 
rity lines, i.e. spin-spin correlators, carry frequency in 
this case. We show that all superconducting parameters, 
including the transition temperature, should depend on 
the properties of the spin system via the spin-spin auto- 
correlation function. In a non-equilibrium state, the spin 
system is described with the help of three correlators (re- 
tarded, advanced and Keldysh). In the quasiequilibrium 
case, these correlators are connected via the fluctuation- 
dissipation theorem. At the end of Sec. II, using analyt- 
ical continuation on the discrete Matsubara frequencies, 
we rederive a relatively simple equation on the supercon- 
ducting transition point in the quasiequilibrium case. 

In Sec. Ill, we address the question of calculating 
spin-spin autocorrelation function in the spin-glass state. 
We propose a method which combines ideas of Thou- 
less, Anderson, and Palmer (TAP)E£l and virial expansion 
method. In the framework of this approximation, the 
A-spin problem is solved exactly while the other spins 
are replaced by a mean value, which plays the role of a 
random mean field. A distribution function for the ran- 



dom mean-field is derived analytically for the case of the 
three-dimensional Heisenberg model with the RKKY in- 
teractions. Let us note that despite the high-temperature 
asymptotics developed in Ref. o, in the glass phase, such 
an expansion does not contain any small parameter. One 
can expect, however, that at large enough N this expan- 
sion gives a quantitatively correct result. At low N, we 
may expect to derive qualitatively acceptable results and 
get some insight into the complicated problem. To il- 
lustrate how the method works, we analytically derive 
the spin-spin autocorrelation function in the first virial 
approximation. 

In Sec. IV, we derive a correction to the supercon- 
ducting quantum critical point. We show that the shift 
of the Abrikosov-Gor'kov's result is quite small (the crit- 
ical concentration of magnetic impurities increases about 
10% compared to the non-interacting case). We also dis- 
cuss the back effect of the superconductivity on the spin 
system. The phase diagram for the system under consid- 
eration is given. 

In Sec. V, we study a non-equilibrium or ageing dy- 
namics in the superconducting transition point. We pro- 
pose an experiment which should provide an explicit 
manifestation of the ageing dynamics in the supercon- 
ducting quantum critical point. If an external magnetic 
field is switched off below T g , the magnetization does 
not disappear immediately. Spins magnetize the elec- 
trons and this effect leads to a stronger suppression of 
superconductivity compared to the case of arbitrarily 
oriented spins. The remanent magnetization logarith- 
mically slowly decreases and drives the system towards 
superconductivity. 



II. NONEQUILIBRIUM GOR'KOV EQUATIONS 

A. The model 

The starting point for the problem is the following 
Abrikosov- Gor 'kov 's Hamilt onian : 



Wag — Wbcs + ^eS, 
where the first term is the usual BCS Hamiltonian: 
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(1) 



WbCS = J |V4 (jj^ - Ipa ~ A^^V^a j <P* 



(2) 



and the second one is the exchange interaction: 



n 



eS 



J Ul £ u (r - r„) (S a a aP ) ^ j d 3 r. (3) 



We neglect finite-size effects and consider the following 
form of the exchange potential: 



u(r) = u a 5 (r) 



(4) 
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The effective spin-spin interactions are described by the 
following Hamiltonian: 



W ss = i^J(r a -r 6 )S a S(, 



In Eq.(§), J(r) is the RKKY interaction: 

cos (2p F r) 
J W = J o ~3 ■ 

The amplitude of this interaction reads: 

mp F 2 



(5) 



(6) 



(7) 



This quantity is connected with the spin-flip scattering 
time as follows: 



J n s T s = [4tt 2 S(S+1)] 1 



(8) 



Let us emphasize that the number in the right-hand side 
of Eq.(@) is very small. 



B. Keldysh Green functions 

In order to find the superconducting transition point, 
we should first derive Gor'kov's equation for the Green 
functions. To treat a possible non-equilibrium dynamics, 
we use the Keldysh technique. Below we recall the basic 
definitions and notations: 

The electron Green functions are defined as follows: 



G>(1,2) 



and 



G<(1,2) = 



V T (1)V>{(2) ^t(1)^(2) 
-^(1)^(2) -^(1)^1(2) 



V}(2)^t(1) V>i(2)V T (i) 
-V{(2)^(1) -V>l(2)V>l(l) 



(9) 



(10) 



Retarded, advanced and Keldysh Green functions are 
then constructed in the following way: 

G R {l,2)=6(t 1 -t 2 ) [g > (1,2)-G < (1,2)1 , (11) 



G A (l,2) = -9(t 2 -t 1 ) G > (1,2)-G < (1,2) , (12) 



and 



G K (1,2) = G > (1,2) + G < (1,2). 



(13) 



It is convenient to introduce a compact notation for 
the 4x4 matrix, (here we use the Larkin-Ovchinnikov 
representational) : 



G(l 



2)= (6*{1,2) G*(l,2)\ 



The Green function satisfies the following matrix equa- 
tion: 



[G.7 1 - E] G = 1. 



(15) 



where Go is the Green function of a normal metal without 
impurities and £ is the self-energy. 
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^ — ^ y G(t; t2 ) 

FIG. 1. Electron self-energy in the leading approximation 
with respect to the magnetic impurity concentration. 7's are 
the vertex matrices [see Eq,(|2l"])]. 

In the case under consideration, we have to find the 
self-energy associated with the scattering on magnetic 
impurities. In the leading approximation on the impu- 
rity concentration, the self-energy has the form shown on 
Fig. 1, where 7's are the vertex matrices (see below). In 
this section, we will consider the case of zero-field cool- 
ing, so that there is no remanent magnetization in the 
spin system and no external magnetic field. In this case, 
we can average Green functions out over the spin degree 
of freedom. In the absence of the RKKY interaction, 
the electron scattering on magnetic impurities is elastic 
in the Born approximation. If we switch on the inter- 
actions between spins, self-energy E becomes non-elastic 
and in some sense analogous to the self-energy due to the 
electron scattering off the phonons. As a result we obtain 
the following expression for the self-energy: 
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±t s S{S ■ 



-^l^a {Tz9fz) ap C\nrfu, 



(16) 



where the Greek indexes label matrix elements in the 
Keldysh space. The structure of the self-energy in 
the Nambu-Gor'kov space is determined by the product 
(f z gf z ) with 



1 

-1 



(17) 



In Eq. ((TJ), we have introduced quasiclassical Green 
function: 



d£ P G(p;ii,i 2 ), £ P = v F (p-p F ) , (18) 



and the following spin-spin correlators: 

C > (h,t 2 ) = C < (t 2 ,h) = -i(S(ti)S(t 2 )) • (19) 
The corresponding Keldysh matrix is defined as 



C 



c k c r 
C A ' ' 



(20) 
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with C R , C A , and C K denned via Eqs.flfi]-|l|). In- 
deed, the spin-spin correlators are proportional to the 
unit matrix in the Nambu space. The vertex matrices 
has the following, explicit form in the Larkin-Ovchinnikov 
representation:^ 







1 
1 



and 7 



1 

1 



(21) 



C. Gor'kov equations 

Since we are not interested in the momentum depen- 
dence of the Green function, it is convenient to integrate 
out the corresponding redundant degree of freedom. Us- 
ing Eq. (|l5|) , we obtain: 

/ d£ p {[G o 1 -t]G-G[Gv 1 -£]}=0. 

As a result, we obtain the following equation which in- 
volves quasiclassical Green function only:E3 

- dg(h,t 2 ) dg(t 1 ,t 2 ) „ 

T~ T- 1 7T. T, 



dti dt 2 

- i[A(ti)g(ti,t 2 ) - g(fa,fa)A(fa)] 

- i [E og -got] (tufa). 



(22) 



where symbol "o" means convolution with respect to the 
time variable. 

Note that each element of g is a 2x2 matrix in the 
Nambu space: 



g 



a fa 
-PI -a 



(23) 



where a is the normal electron Green function and j3 is 
the Gor'kov (anomalous) function. In the absence of the 
external currents and/or magnetic field, Pi = Pi- Note 
also that matrix ( |l4| ) is a matrix in the direct product of 
independent Keldysh (time-reversal) and Nambu spaces. 
It is a question of convention how to construct this ma- 
trix. One can use either form (jlj) or another way: 



a P 
-p* -a 



(24) 



with d and P being matrices in the Keldysh space. Since 
we are looking for the transition point, we should first de- 
rive equation linear on the superconducting Green func- 
tion p. Thus, form ( g4j) is more convenient for our pur- 
poses. Using the notations introduced above, the order 
parameter A can be written as 



A(t) = ir\\\vp K (t,t + 0), 



(25) 



where v is the density of states per spin at the Fermi- 
surface. In the matrix notations, the order parameter 
takes on the form: 



A(t) 



A(t) 
-A*(t) 



(26) 



In the Keldysh space, A is proportional to the unit ma- 
trix (we neglect superconducting fluctuations). 

To find the superconducting transition point we have 
first to extract the equation on the Gor'kov Green func- 
tion p. After some algebra, we obtain: 



_d d_ 

dTi~dT 2 

i 



At s S{S 



P^ v (ti,fa) + i[A(ti) + A(fa)] (tufa) 

ryj {l»olp6 [ a <y P C\ v o p Sv + P„ p C\ n o asu) (h,t 2 ) + [a^s ° C\ v p ap + PpS ° Cx^p] £2)7^7^} , (27) 



Although, the spin system may be out of the equi- 
librium state, the electron system may be considered in 
equilibrium. This yields: 

g K = g R °.f-fog A , 

where / is the Fermi-distribution [/(e) = tanh (Sf) i n 
the energy representation.] 

Let us now separate slow and fast dynamics in the spin 
correlator by performing Wigner transformation: 



sp R (s) 



C(u,t)= d(ti-t 2 )C(h-t 2l t)e-^- t2 K (28) 



In Eq. (pg)) , parameter t = (ti +fa)/2 describes slow dy- 
namics and, in some sense, is similar to the waiting time. 
From Eqs.(Q|6|,^2^^^|g|^), we obtain the follow- 
ing equations: 
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> - /■ 1 . '//-x d 2 a R (e) 
ACt)a R (e) + l A"(t)—^l 



4t s S(S+1) J 2tt 

— OG 

:[c k (lo) [a R (e)p R (e - u) + P R (e)a R (e - to)] +C r (lu) [a R (s)p K (e - lu) + p R (e)a K (e - lo)]} (29) 



and 
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el3 A (e) + 



A(i)a» + ~A"(i) 



i d 2 a A (e) 
4 ( > de* 



+ 00 



At s S(S+1) 



dco 
2^ 



: {c k (lu) [a A (s)p A (s - u) + (3 A {e)a A {e - u)] + C A (cu) [a A {e)fi K {e - u) + (3 A {e)a K {e - u)] }. 



(30) 



The corresponding Keldysh function reads: 

e 



(3 K (e) = [f3 R {s) ~l3 A {e)] tanh 



2T 



(31) 



Note, that all quantities in Eqs.( |29|j30| ) weakly depend on 
the "waiting time." Together with Eq.(j25j), Eqs.(||Q 
form a closed equation system, which contains all the in- 
formation about the superconducting properties of the 
physical system under consideration. The spin system 
is completely described by the three functions: C R , C A , 
and C K . These functions are not necessarily connected 
with each other in an out-of-cquilibrium state and depend 
upon the details of the state and history. 

At the superconducting transition point, Eqs.(^9 30) 
can be simplified with the aid of the following relations: 



l (e) = -a A (e) = 1 



and. thus, 



a K (e) = 2 tanh 



2T 



(32) 



(33) 



D. Quasiequilibrium regime 

Let us consider a regime when the dynamics of the 
spin glass state has two well separated time scales: The 
first one is "waiting time" i, which is comparable with 
the typical time of an experiment. At these time scales, 
large energy barriers change. The second typical time 
is fi/Tg. At these time scales, a quasiequilibrium state 
is reached within the large energy barriers. Tunneling 
processes through the large barriers are supposed to be 
very rare at such a time-scale. In this approximation, the 
waiting time can be considered as an independent exter- 
nal parameter and equilibrium techniques can be used. In 
this equilibrium regime the fluctuation-dissipation theo- 
rem holds and it significantly simplifies the calculation. 

In the equilibrium, Keldysh, retarded and advanced 
spin-spin correlators are not independent functions any- 
more. They are connected via the fluctuation-dissipation 
theorem, which follows directly from the Gibbs's distri- 
bution and in our language can be formulated as follows: 



C K {u)= [C R (cj)-C A (u;)]coth[^] 



2niS (uS) 



E 



(34) 



|S|j)| 2 e-*, 



where \ k) is a quantum state in the spin system and is 
the corresponding energy level. Note, that for spin corre- 
lators, the last "static" term in Eq.([34j) does not vanish. 

In the equilibrium, it is convenient to introduce Mat- 
subara Green functions f3{e n ) and a{e n ) which depend 
on the fermion frequency e n = ir(2n + 1)T and the Mat- 
subara spin-spin correlator: 



l/T 

C (w m ) = J dT (T T [S(t)S(0)]> e~^ T , 



(35) 



where 



S(r) = e- HT S(0)e Hr 
and uj m = 2-KmT is the bosonic Matsubara frequency. 
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FIG. 2. Contour in the complex plane to used to perform 
analytical continuation on discrete Matsubara frequencies for 

En > 0. 

Using the contour shown on Fig. 2, one can per- 
form analytical continuation on the discrete Matsubara 
frequencies.EJ To do this we first note that Keldysh spin- 
spin correlator consists of two different parts: one is pro- 
portional to 5(ui) and the other, we denote it as C K , sat- 
isfies the relation C K '(«) = [C r {uj) - C a {uj)] coth [^] . 
The former part gives an elastic contribution to the self- 
energy. It is very easy to treat this term, we just replace 
8{uS) by a delta symbol y<L m o and the integral ( p9| ) by 
the sum over u> m . The latter part is more complicated. 
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As an example, let us consider the first and third terms 
in the right-hand side of Eq.(^). After analytical con- 
tinuation on discrete Matsubara frequencies e„ > 0, we 
get for these terms: 




C r (w)-C a {uj) 



+C r (lo) [(3 R (ie n -Lu)-(3 A (ie 
dm ~ 



;)] tanh 



coth 



2T 



=-C R (cu)p A (is n 
Zn 



^C A (u)0 R (is r , 

Z7T 



2T 

- uj) coth 

- uj) coth 
uj) coth 



we propose a method of calculating this correlator which 
introduces the concept of a random mean field. This idea 
is similar, to the one used by Thouless, Anderson and 
Palmer £3 We combine the random mean- field method 
with the virial or cluster expansion [see e.g. Ref. b|. 

Let us consider a spin in the mean-field of the other 
spins h = ^Jab^a- The distribution function for this 

a 

random magnetic field can be naturally defined as fol- 
lows: 



2T 



2T 



P[h] = 

or, equivalently, 
d 3 \ 




(38) 



j.s 



2T 



P[h] 



(36) 



(27T) 



exp 



-ah 



Jab^a 



J.S 



r 3 



We have used Eqs.( p2|3^ ) valid at the transition point 
and also the relation tanh [ "^" 1 = — coth [^] . It is 

easy to see, that in Eq.(|3^), all C and f3 functions in- 
volved are analytical inside the contours of integration. 
Thus, the total integral is determined by the poles of 
the cotangent function only. As a result, the whole ex- 
pression ( |36| ) can be rewritten as a sum over the bosonic 
Matsubara frequencies: 







Similar procedure can be done with the second and 
forth terms in Eq. (p9|) . Summarizing, we recover the 
quasiequilibrium equation on the superconducting tran- 
sition point. 



\f3(s r , 



T 



2t s S(S+1) 

x {(3(e n - uj m ) + /?(£„) sgne„sgn ( 



-^C(u; m ) 

■Wm)} (37) 



Note that the equilibrium expression is the same in para- 
magnetic and spin-glass phases. However, the spin-spin 
autocorrelation functions are different in these phases. 



III. SPIN-SPIN AUTOCORRELATION 
FUNCTION 

A. Random mean-field 

To find the superconducting transition point, we have 
to solve Eqs. (Hg) or Eq.@, which requires knowing 
the spin-spin autocorrelation function(s). Calculation of 
this correlator for the quantum Hcisenberg spin-glass is a 
very hard problem even in the quasiequilibrium case. To 
get some insight into the problem under consideration, 



(39) 



where the averaging over the random spin-spin interac- 
tion and over spin orientations is implied. For the RKKY 
interaction (^), averaging means: 



</(</)> 



f 

V 



d 6 r 



d<f> I Jq cot < 



2tt 



(40) 



where V is the volume of the system and / is an arbitrary 
function. 

From Eqs.(^,E3), we obtain: 



P[h] 



d 3 X 

(27T) 3 



iXh 



n 



4tt Jq 




(41) 



Performing averaging over the spin orientations and eval- 
uating the corresponding elementary integral we get the 
distribution function for the random mean field: 



P[h] = \ 

7T 2 ( a 2 



h 2 Y 



a = —J n s S, 



(42) 



where n s is the concentration of magnetic impurities and 
S is their spin. 

Let us note that for a different model of spin-spin inter- 
actions the distribution function (|3^) would be different. 
For example, if we start with the Gaussian distribution of 
spin-spin couplings in the Shcrrington-Kirpatrick model, 
the distribution function of the mean field is necessarily 
Gaussian as well. 

In the non-equilibrium case, the simple calculation pre- 
sented above is not valid. The distribution function of the 
mean-field should depend on time and on the type of the 
non-equilibrium state. For example if we study a non- 
equilibrium dynamics due to a slowly decaying remanent 
magnetization, one should modify the definition of the 
distribution function as follows: 



G 




P[h, l] = Mh-^ J Qb S a ^J x 6 (m(t) -J2 S * 



where m(£) is the remanent magnetization. 

B. Virial expansion 

Now, let us consider the following set of Hamiltonians: 
Hi = Sh, (44) 



Ti-2 — ^12^182 + Sihi + S2I12, 



k 

n {k) = -J2 Ja fc s a s fc + ^s a h a , 

a, 6 a 



(45) 



(46) 



where both J's and h's are random and distributed ac- 
cording to Eq.(|40|) and Eq.(ff2|), correspondingly. 

In principal, one can calculate any quantity using 
Hamiltoninans (|4|) (one spin) , (|||) (cluster of two spins) , 
etc. and than average out the quantities of interest using 
Eq.([40|) and distribution ((42|). This procedure generates 
a series which can be called virial or cluster expansion. It 
is similar to the virial expansion in the theory of liquids 
and gases. Let us also note, that comparison with the 
Sherrington-Kirpatrick Ising model shows that the first 
virial term is equivalent to the replica symmetric solu- 
tion of the model. After a simple calculation one can 
obtain the equation for the Edwards- Anderson order pa- 
rameter in the replica symmetric case. It is well known 
that the corresponding solution is not stable. The usual 
practice is rto apply the mechanism of replica symmetry 
breaking.Lrtj However, replica technique hardly can lead 
to some explicit results especially in the dynamical prob- 
lems. Moreover, it involves a procedure of analytical con- 
tinuation on replica indexes. This procedure is usually 
poorly justified, since, the behavior of functions to be an- 
alytically continued is typically not known at large replica 
indexes. It is not clear whether the virial expansion can 
help to solve the underlying problems. Technically it 
is quite straightforward and physically very transparent. 
Even though, the expansion does not contain any small 
parameter at low temperatures, one can try to calculate 
a quantity of interest up to some large N (number of 
spins in the cluster) in order to sum up the series. Even 
if the series is not convergent, one can use Pade-Borel 
technique to do the summation. 

If we are interested in studying a system in a non- 
equilibrium regime, the Hamiltonian language described 
above is not appropriate. One has to reformulate the 
problem using the corresponding action on the Keldysh 
closed time contour and do the virial expansion within 
this formalism. This idea will be developed elsewhere. 



C. Matsubara spin-spin correlator. Equilibrium case. 

To give an example of how this method works, let us 
calculate the spin-spin autocorrelation function in the 
first virial approximation. We will consider quasiequi- 
librium regime only. For the Matsubara correlator, we 
obtain: 



C (w m ,/i) 



E K fc il s (°)l fc 



J/l 2 x 



k\ ,k 2 



e k 2 - (k t 



[ e -^ 



(47) 



For spin S = 1/2, in the first virial approximation, we 
get: 



C(uj m ,h) = ^<5 Wm , 



h 



h 2 



tanh 



2T 



(48) 



Note, that for an arbitrary spin value S, the tangent fac- 
tor is replaced by the correspon ding Brillouin function. 
We have to average out quantity ( f48|) with respect to the 
distribution function (B2): 



C(w m ) = / d 3 hP[h]C(u m ,h). 



Evaluating the corresponding integral, we obtain 
C(uj m ) = —S Umt0 + SC(uj m ), 



(49) 



(50) 



where 



SC(u> m ) = 4ai 



2a 

7T 



res ~ res 

z—ia z—i\uj 17l \ 



I(a,uj m ) + 2a 



dl(a,w m ) 



with 



9{*) = 



(z 2 + ui 2 ) (z 2 + a 2 Y 



da 



■ tanh 



(51) 



2T 



and 

I(a,w Tl 



: Re 



1 1 . iu m \ 



ia \ 
2rfT J 



At a low temperature T -C a ~ T g , Eq. ( plj) can be simpli- 
fied and we obtain the following zero-temperature result: 



SC(uj) 



2a a 2 ~aj 2 +uj 2 In 

7T {L0 2 ~a 2 f 



(52) 
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Let us emphasize the following property of correlators 
© and (H): 



lim C(aj m ) = S(S+ 1) <L,„o, 



(53) 



which can be verified by a straightforward calculation. 
This property means that if we switch off the interac- 
tions between the spins, the corresponding autocorrela- 
tion function turns into the time-independent correlation 
function of a free spin. 

In the second virial approximation we should calculate 
correlator (|47| ) using eigenfunctions of Hamiltonian (45). 
There are four quantum states possible in this case. Even 
at zero temperature, finding the ground state energy re- 
quires solving spectral problem which turns out to be an 
equation of the fourth order. Moreover we should average 
out correlator ( |47| ) over the random fields and coupling 
constant: 



C(cu m )= / d 8 hid 3 h a P[h 1 ]P[h2]<C(t«; TO ;J s hi,h2)) J 



(54) 



This is hard to handle this analytically. Numerical work 
is required in this case. 



IV. SUPERCONDUCTING TRANSITION POINT 

A. Quantum critical point. Equilibrium case. 

To find the superconducting transition point, we have 
to solve Eqs.(|2^, ^) or Eq.(|37j) in the quasiequilibrium 
regime. Let us first consider the latter case. At a fi- 
nite temperature, Eq.(|37|) is an infinite system of coupled 
equations corresponding to different Matsubara frequen- 
cies. As temperature goes to zero, the picture is simpli- 
fied and Eq.(p7|) turns into an integral equation: 



N + - 



(3(e) - A = - 



i r dw 

2t s S(S + 1) J 2^ 



x - s(s + i)<y(w)] 

x \J3(e + uj) + (3(e)sgne sgn (e + u))] . (55) 

In the absence of spin-spin interactions the anomalous 
Green function (3 reads: 



/3 (0) (e) = 



A 



(56) 



Using definition (p5[), multiplying Eq.([55|) on (}( >(e), and 
integrating it over e, we get: 



■kv | A | 



de 
2^ 

du> 
2^ 



2t s S(S+1) 
P(w) -S(S + 1) 5{uj)] (57) 



de 
2ttA 



u) + (3(e)sgne sgn (u + e)} /5 (0) (e). 



As it will be shown below, this equation can be solved 
using the perturbation theory. Due to numerical reasons, 
the term in the right-hand side can be treated as a small 
perturbation. Thus, one can use expression ( |5^ ) for the 
Green functions in Eq.([57j). In this case, the integral over 
e in the right-hand side of ( |57j ) can be easily evaluated 
and we get the following expression for the integral: 



2t s I{ujt s ) = 2t s 



1 + \co\t s 



\lu\t s (1 + Mt,/2) 



ln(l+|w|T,U (58) 



Correlator C follows from Hamiltonian (||) and distri- 
bution function defined by relation (phf). There is only 
one dimensional parameter upon which it can depend: 
a ~ Jo n s- Thus, in the equilibrium, the correlator is 
bound to have the following form: 



CM 



a V a 



This statement is true both in the high-temperature 
and low-temperature regions. Let us make the follow- 
ing change of variables in the integral over to in Eq. (^7|) : 
u> — > x = uj/a. Then, this equation can be rewritten as 
follows: 



In 



2 1 



AnS(S+l) 



dx 



C{x)-S{S+l)8{x) I(ax), (59) 



where function I(x) is defined in ( p8[ ) and a — ar s . 
The latter quantity is just a number. From, Eq.(^|) and 
Eq.©, we find 



1 



6?r(5+i; 



< 1. 



One can see that this number is very small. If, for- 
mally, we put (7 = 0, we will get zero in the right-hand 
side of Eq.(|9|) [see also (|||)] and recover the Abrikosov- 
Gor'kov's formula for the quantum critical point: 



27 

= — » 1.13. 
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If we keep a finite, the right-hand side term in Eq.(|5£ 
will give a correction to the quantum critical point due to 
the spin-spin interactions. For example, we can use corre- 
lator ([32]) calculated in the previous section. Evaluating 
the corresponding integral numerically we get correction 
to the critical concentration of magnetic impurities con- 
centration: 



Sn 



(i) 



,(AG) 



0.073. 



'.sO 



Similarly, we can obtain the shift corresponding to the 
second virial approximation (p4|). This has been done 
numerically. The corresponding result is 



,( 2 ) 



JAG) 



0.085. 



'.sO 



We see that these corrections are quite small. This small- 
ness is connected with the parameter T g r s present in the 
theory and should exist in any order of the virial expan- 
sion. Thus, we conclude that the shift of the Abrikosov- 
Gor'kov quantum critical point is about 10%. 



B. Equilibrium phase diagram 

Another interesting question which is worth discussing 
is the effect of the superconductivity on the effective spin- 
spin interactions. jlhis questions has been addressed long 
ago \w- Abrikosovila Recently it was revisited by Aristov 
et al.t3 The effective spin-spin coupling is determined by 
the following expression: 



J(r) = 2u§T£[|a(e n ,r)| 



>r)r 



(60) 



Tco 



Normal Metal 




High-T expansion 



Superconductor 



Low-T expansion 



Paramagnet 
Spin glass 

f-3 n s o 

FIG. 3. Schematic phase diagram for a superconductor 
with interacting magnetic impurities. The dashed curve is 
the Abrikosov-Gor'kov's transition line. 

In the absence of superconductivity, we expect the spin 
glass state to survive even at very low concentrations of 
magnetic impurities if T — 0. In a superconductor, the 
spin glass state should disappear at 



r 



T C Q 



We see that not only spin-spin interactions change the su- 
perconducting phase diagram, superconductivity in turn 
suppresses magnetic ordering at low temperatures and 
spin concentrations. 

Combining these results with the ones obtained in the 
previous sections, we can plot the phase diagram for the 
physical system under consideration (see Fig., 3). Let 
us note that the high-temperature asymptoticsB matches 
the low-temperature one at a temperature T* ps 2Jqu s , 
which is approximately the temperature of the spin-glass 
transition. 



where a(e„,r) and /3(e„,r) are Green functions in the 
real space. A straightforward calculation yields the fol- 
lowing form of the effective RKKY interaction below the 
superconducting transition point: 



J(r) = 



cos (2kpr) 



^2 ^ AF 



-2r/£ 



(61) 



where £ = vp /A and antiferromagnetic contribution 
Jaf oc A/cf is quite small. 

At the transition line, £ = 00 and A = 0, so we re- 
cover the usual RKKY formula. However, well below the 
transition the interaction becomes screened. This fact 
has a very transparent physical explanation. The RKKY 
interactions appears due to the fact that a spin polarizes 
the normal electrons around it. The Friedel's oscillations 
in the electron density lead to the Ruderman-Kittel os- 
cillations in the spin-spin coupling. In a superconductor, 
the spins are not polarized in the ground state. Thus at 
the distances larger than Cooper pair size £, the indirect 
spin-spin coupling must be suppressed. 



V. AGEING IN THE SUPERCONDUCTING 
QUANTUM CRITICAL POINT 

Let us study the following experiment. Consider a su- 
perconductor at a low temperature with the concentra- 
tion of magnetic impurities n s < n s0l where n„o is the 
equilibrium critical concentration found in the previous 
section. Now, let us switch on a large magnetic field. 
This field polarizes the magnetic impurities and elec- 
trons and destroys superconductivity. Then, we switch 
off the magnetic field. The electron system equilibrate 
very quickly. However, in the system of magnetic impu- 
rities, one finds a remanent magnetization which decays 
very slowly with time. This magnetization acts on the 
electrons as an external magnetic field, which affects only 
the spin degree of freedom but not the orbital one. In 
some sense the effect due to the remanent magnetiza- 
tion on the electron system is equivalent to the effect 
of an external magnetic field applied parallel to a two- 
dimensional superconducting sample. This leads to the 







further suppression of superconductivity. Thus, even af- 
ter the external magnetic field is switched off, the super- 
conductivity may be absent in such an experiment due to 
the ageing effects in the system of magnetic impurities. 

To get some qualitative insight into the non- 
equilibrium case, let us consider a small enough remanent 
magnetization and neglect corrections to the Abrikosov- 
Gor'kov form of the spin-spin correlator. In this case, 
we obtain the following set of equations which describes 
the superconducting transition and determines the order 
parameter: 

fe - i6/A /3(e) = Aafe) - — a(e)/3(e), (62) 

where a(e) and (3{e) are superconducting Green functions 
subject to the following constraint: 



a 2 (e) + \/3(e)\< 



1 



(63) 



and the order parameter is defined as usual: 

A = 7rHA|T^/3(£„). (64) 



In Eq.(|62j), we have introduced the splitting of the Fermi- 
surface due to the remanent magnetization: 



Sfi = u m(t), 



(65) 



where uq is the integral over the exchange potential [see 
Eq.(||)] and m(t) = n s (S) is the remanent magnetization 
which slowly depends on time. Let us note, that upon 
m(t), all superconducting quantities in Eq.(|62|) acquire 
similar slow timc-dependence. 

From Eqs.(62 — [34]), it follows that for large enough 
values of remanent magnetization (5/i ~ r s _1 , the super- 
conducting transition becomes of the first order. The 
order parameter jumps from the zero value up to a finite 
value. Let us consider smaller values of the magnetization 
and find the corresponding second-order superconducting 
transition point. From Eqs.(^| — |64|), we get: 



1 =TTV\\\T C J2 



\e\ - i (S/j,) sgne + r s 



(66) 



At zero temperature, the corresponding integral can be 
easily evaluated and we get: 



1 = In 



27 y/SjT 2 



Finally, we obtain the following time-dependent critical 
concentration: 



n s (t) = ^Jn 2 s0 - cm 2 (t), 
where constant c has the form: 



(67) 



S (S + 1) mp-p 



(68) 



The dynamics of the remanent magnetization in the spin- 
glass state can be usually well-fitted via the following 
formulary 

m(t) = mo — vhit. 

Where mo is of the order of the remanent magnetization 
just after the magnetic field is turned off and quantity 
v is called magnetic viscosity. Thus, the superconduct- 
ing critical point is slowly flowing towards its equilibrium 
value: 



n s (t) = n s0 



- [mo - v\iit] z 



(69) 



We see that after some time, which can be macroscopi- 
cally large, the superconductivity should appear again. 

Let us note that non-equilibrium effects at the super- 
conducting transition point does not necessarily have to 
be connected with the dynamics of remanent magnetiza- 
tion. Another possible experiment could be done as fol- 
lows. Let us consider a superconductor with the concen- 
tration of magnetic impurities such as n]^ G ^ < n s < n s o, 
where again n^ G ^ is the Abrikosov-Gor'kov's critical 
concentration at zero temperature and n s o is the real 
critical concentration at zero temperature with the ac- 
count of interactions. Consider large enough initial tem- 
perature so that no superconductivity is present at the 
beginning. Then, let us cool a sample down very quickly 
at zero field. One should expect that superconductiv- 
ity appears only some time after the sample was cooled 
down due to the slow relaxation processes in the spin 
subsystem. 



VI. CONCLUSION 

We have considered a superconductor with interacting 
magnetic impurities. The central result of this paper is 
the dependence of the superconducting properties in such 
a system upon the spin-spin autocorrelation function. At 
low temperatures, one can expect to observe ageing ef- 
fects in the superconducting transition point. The limit- 
ing manifestation of the ageing would be the observation 
of the spontaneous appearance of superconductivity after 
some macroscopic waiting time due to a slow change in 
the remanent magnetization in the spin subsystem. 

In the equilibrium, we predict that the Abrikosov- 
Gor'kov's critical line shifts towards higher temperatures 
and impurity concentrations due to the spin-spin interac- 
tions. Let us note that the high-temperature expansion 
diverges as T — > and, formally, predicts very large shift 
of the critical point. The two asymptotics match at a 
temperature T* « 2Jon s . This gives a good estimate 
for the temperature of the spin-glass transition T g . At 
T = T g , a crossover from the high-temperature (param- 
agnetic) asymptotics to the low-temperature (spin-glass) 
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one takes place. We have shown that the actual shift of 
the superconducting transition due to the RKKY inter- 
actions is small at any temperatures. This is connected 
with the small parameter T g r s which exists in the theory. 

Let us also mention that our calculations are valid only 
for temperatures T 3> Tk, where Tk is the Kondo tem- 
perature. In this case we can neglect higher order pro- 
cesses in the electron scattering off the impurities. At 
T ~ Tk the magnetic impurities become screened, which 
presumably should lead to the further enhancement of 
superconductivity. 

Let us note that some deviations from the Abrikosov- 
Gor'kov's curve have been observed experimentally.Q The 
measurements showed some increase in the critical con- 
centration at low temperatures. Except the effects dis- 
cussed in the present paper, this increase can also be con- 
nected with the phenomena similar to the ones discussed 
in Refs. plj^ Namely, optimal fluctuations in the distri- 
bution of magnetic impurities may lead to the formation 
of superconducting islands coupled via the Josephson ef- 
fect. This effect may also lead to some increase in the 
critical concentration. However, we do not expect this 
increase to be very large due to the suppression of the 
Josephson coupling by quantum fluctuations £3 To iden- 
tify which effect is dominant at low temperatures, addi- 
tional experiments are highly desirable. In either case, we 
expect that some interesting non-equilibrium phenomena 
(hysteresis, ageing, etc.) should reveal themselves at low 
temperatures. 
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